COHOMOLOGY OF CONGRUENCE SUBGROUPS OF SL 4 (Z) 
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Abstract. Let N > 1 be an integer, and let L = Tq(N) C SL^Z) be the sub- 
group of matrices with bottom row congruent to (0, 0, 0, *) mod N. We compute 
H 5 (T;C) for a range of N, and compute the action of some Hecke operators on 
many of these groups. We relate the classes we find to classes coming from the 
boundary of the Borel-Serre compactification, to Eisenstein series, and to classical 
holomorphic modular forms of weights 2 and 4. 



1. Introduction 

1.1. Let n > 1, and let T be a congruence subgroup of SL n (Z) of level N . Let Sn 
be the subsemigroup of the integral matrices in GL n (Q) such that (T, Sn) is a Hecke 
pair. 

We denote by Jf?(N) the C-algebra of double cosets TSnT. This algebra acts on 
the cohomology and homology of T with any coefficient CSV-module. When a double 
coset is acting, we call the map defined by its action a Hecke operator. In this paper 
we will work only with the trivial coefficient module C 

Let I be a prime not dividing N, and let D(l,k) be the diagonal matrix with 
diagonal (1, ... ,1,1, ... ,1), where the number of Ts is k. Then J4?(N) contains all 
double cosets of the form TD(l, k)T, and we denote the corresponding Hecke operator 
by T(l,k). Fix a prime p not dividing N, and an embedding of Q p into C. Let 
Go = Gal 



Definition 1.2. Let V be an J^(iV)-module and suppose (3 G V is an eigenclass 
for the action of J$?(N). For I prime to N, write T(l,k)(f3) = a(l,k)(3 , where the 
a(l,k) G C, k = 0,...,n are algebraic integers. Let p be a continuous semisimple 
representation p : Gq — > GL„(Q P ), unramified outside pN, such that 



.l)*/fc(*-i)/2 a (/ j k ) X k = det(/ - p(Frob i )X) 
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for all I not dividing pN. Then we shall say that p is attached to (3. 

For example, theorems of Eichler, Shimura and Deligne imply that if n = 2 and 
V is the Hecke- module of classical holomorphic modular cuspforms for SL 2 (Z), then 
there always exists a p attached to any Hecke eigenform /3 G V. 

Standard conjectures (for example in |T5|]) state that if V is the cuspidal cohomology 
of T with trivial coefficients, then any Hecke eigenclass should have an attached p- 
adic Galois representation. In fact, one has the stronger conjecture that there should 
be an attached motive (cf. We can also extend the conjecture of [|l^] to include 



all the cohomology — in principle, the theory of Eisenstein series should allow one to 
reduce the extended conjecture to the one for cuspidal cohomology. 

In a series of papers || §, [l], |1X| , the first author with a number of coworkers has 
tested this conjecture and a mod p variant of it when n = 3. Other tests for n = 3 can 
be found in the work of van Geemen and Top with van der Kallen and Verberkmoes 



26, 27, 28]. The purpose of this paper is to make the first computational tests of 



this conjecture for n — 4. 



1.3. We work here with T = Tq(N), defined as the subgroup of SL n (Z) consisting 
of the matrices with last row congruent to (0, ... , 0, *) modulo N. 

Conceptually, there is no big difference in computing the cohomology of T between 
the cases of n = 3 — or even n = 2 — and the case of general n. One needs to write 
down a simplicial complex C homotopic to the chains of a _Br(l)-space, and compute 
the cohomology. 

When n = 2 or 3, the most interesting part of cohomology, namely the cuspidal 
part, occurs in the top dimension of C. Therefore to compute it all we have to 
compute is the cokernel of a coboundary map. However, when n = 4 the cuspidal 
cohomology occurs in dimensions 4 and 5, whereas the virtual cohomological dimen- 
sion of T, and hence smallest possible dimension of C, is 6. So now the cohomology is 
a subquotient, which adds considerably to the complexity of the computer programs. 
The cuspidal cohomology in degree 4 is dual to that in degree 5, so we concentrate 
on computing the latter in this paper. 

But the big difference between n = 2, 3 and n = 4 occurs when we try to compute 
the Hecke action on the cohomology. In the top dimension we can use the Ash- 
Rudolph algorithm or its variants ||, as was done for n = 3 in the works cited above. 
However, for n = 4, where we look just below the top dimension, a brand new idea 
was necessary. This is due to the second author, and is the subject of Thus 



this paper is also a test of the algorithms proposed in Fj| , and they pass with flying 
colors. It is an open question whether the algorithms of [0 terminate in a finite 
number of steps; in practice, though, they always terminate quickly, and we used 
them here without problems. 
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Our computations of H 5 {T (N); C), detailed in §|5|, were made for N < 53.Q No 
cusp forms were discovered, but some interesting phenomenology of the boundary 
cohomology was observed, as discussed in §|B|. This leads to some open questions 
about the cohomology of the boundary, which are discussed there. Since we have not 
completed the Hecke computations for some high levels near 53, it is possible that we 
do have a cusp form that we haven't yet identified as such. There is a bound due to 



Fermigier [T6] that states that if N < 31, there cannot be any cuspidal cohomology. 

For high levels our computations produced large sparse matrices, as large as 110464x 
30836 for level 48. To perform linear algebra with these matrices, we used a version 
of the Lanczos algorithm mod p, in the spirit of LaMacchia-Odlyzko [[0J (see also 



1.4. We thank Eric Conrad for some assistance with programming, and Peter Woit 
for excellent computing support. We thank David Ginzburg for conversations at the 
beginning of this project. 



2. Background 

2.1. Let V be the 1R- vector space of all symmetric n x n matrices, and let C C V be 
the cone of positive-definite matrices. Then the group G = SL n (M) acts on C on the 
left by (g, c) \— > g ■ c • g l , and the stabilizer of any given point is isomorphic to SO n . 

Let X be C mod homotheties. The G-action on C commutes with the homotheties 
and induces a transitive G-action on X. The stabilizer of any given point of X 
is again isomorphic to SO n . After choosing a basepoint, we may identify X with 
the global riemannian symmetric space SL n (R)/SO n , a contractible, noncompact, 
smooth manifold of real dimension d = n(n + l)/2 — 1. 

The group SL n (Z) acts on X via the G-action, and does so properly discontinu- 
ously. Hence if T C SL n (Z) is any finite- index subgroup, the quotient T\X is a real 
noncompact manifold except for at most finitely many quotient singularities. We may 
then identify the complex group cohomology H*(T;C) with H*(T\X;C). Although 
the dimension of T\X is d, it can be shown that H t (T\X] C) = 0ifz><i — n + 1 
|T2| , Theorem 11.4.4]. The number v = d — n + lis called the virtual cohomological 
dimension of T. 

In this paper we will always take T to be the congruence subgroup To (A) of matrices 
whose last row is congruent to (0, ... , 0, *) mod N. 

2.2. Recall that a point in Z n is said to be primitive if the greatest common divisor 
of its coordinates is 1. In particular, a primitive point is nonzero. Let 8? C Z n be 
the set of primitive points. Any v e written as a column vector, determines a 



Actually, to avoid numerical instability in floating-point computations, we replace C with the 
finite field F31991 (cf. §£J). 
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rank-one symmetric matrix q(v) G C by q(v) = v ■ v l . The Voronoi polyhedron IT is 
the closed convex hull of the points q(v), as v ranges over 

Note that, by construction, SL n (Z) acts on II. The cones over the faces of II form 
a fan "V that induces a T-admissible decomposition of C [0, p. 117]. Essentially, this 
means that V acts on Y] that each cone is spanned by a finite collection of points 
q(v) where v G and that there are only finitely many T-orbits in "¥ . The fan V 
provides a reduction theory for C in the following sense: any point x G C is contained 
in a unique a(x) G Y, and the set {7 G SL n (Z) | 7 • a(x) = &(x)} is finite. 

2.3. We summarize facts about the well-rounded retract of || |J. There is a defor- 
mation retraction C — > C that is equivariant under the actions of both T = SL n (Z) 
and the homotheties. Its image modulo homotheties is the well-rounded retract W 
in X. The well-rounded retract is contractible, since it is a deformation retract of the 
contractible space X. Hence the cohomology of T\X with coefficients in C is canon- 
ically isomorphic to the equivariant cohomology iff (W; C) where T acts trivially on 
the coefficient module C. This is in turn canonically isomorphic to the complex co- 
homology H l (T\W] C), since C has characteristic zero, which moreover is isomorphic 
to H l (T;C). We will focus on computation of the equivariant cohomology. The di- 
mension of W equals the virtual cohomological dimension u, and the quotient T\W 
is compact. 

The well-rounded retract W is naturally a locally finite cell complex, the cells 
being convex polytopes in V. The group SL n (Z) preserves the cell structure, and 
the stabilizer of each cell in SL n (Z) is finite. The theory of cores and co-cores in 
]7], Chapter 2] shows that the cells in W are in a one-to-one, inclusion-reversing 
correspondence with the cones in the Voronoi fan "V . By abuse of notation, the cell 
in W corresponding to a will still be denoted a. 

2.4. One can give a more precise description of the combinatorics of the cones in Y 
and the cells in W . To each a G Y , we define the set 

M(a) = {v G & I q(v) is a vertex of the face of II generating a.} 

We associate the same set M(a) to the corresponding cell in W, and call M(a) the 
set of minimal vectors of a (because of how W is constructed in 0, [|). Since IT is 
the convex hull of the q(v)'s, it is clear that a — > M(a) is an inclusion-preserving 
(respectively, inclusion-reversing) bijection between the cones in "V (resp., cells in W) 
and a collection of finite subsets of . In principle, this reduces the study of the 
combinatorics of Y and W to the study of the image of /1. For instance, face relations 
r C a in "V are read off from subset relations M(r) C M(a). In practice, determining 
the image of \x calls for explicit computations with real quadratic forms, computations 
whose difficulty grows exponentially as a function of n. The computations have been 
carried out completely for n < 5 by various authors. 
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For the rest of this subsection, we set n = 4 and give more details. We state 
results for the well-rounded retract; these imply their analogues for "f . The image 
of n was computed independently by an d (in essence) [EH]. The cells of W fall 
into eighteen equivalence classes modulo SL n (Z). Let T (for "type") be a variable 
running through these eighteen classes. This partitions the set of cells of W into 
eighteen pieces called the Wt- Any a G Wt is said to be of type T . In each Wt, we 
fix one representative cell <7t, the standard cell of type T. The M(ctt)'s are written 
down explicitly in [21|.Q This determines the image of //, since the image is the union 



of the SL n (Z)-translates of the eighteen M(<tt)'s. 

2.5. To compute the action of the Hecke operators on cohomology, the well-rounded 
retract is insufficient, since the operators do not act cellularly. To ameliorate this, 
we use the sharbly complex. The material in this subsection closely follows 0. 

Definition 2.6. || The sharbly complex is the chain complex {S*, d} given by the 
following data: 

1. For k > 0, S k is the module of formal Z-linear combinations of basis elements 
u = [v i, . . . , v n+k ], where each t>j G mod the relations: 

(a) If r is a permutation on {n + k) letters, then 

[ui,... ,v n+k ] = sgn(r)[r(fi),... ,r(v n+k )\, 

where sgn(r) is the sign of r. 

(b) If q = ±1, then 

[qyt, v 2 ... , v n+k ] = [v u ... , v n+k }. 

(c) If the rank of the matrix (v\, . . . , v n+k ) is less than n, then u = 0. 

2. The boundary map d: S k — > S k -i is 

n+k 

[«!,... , W„+fc] I ► ^(-l) l [wi, ... , Vt, . . . , Vn+k]. 

i=l 

The basis elements u = [f i, . . . , f„+fc] are called k-sharblies. By abuse of notation, 
we will often use the same symbol u to denote a fc-sharbly and the fc-sharbly chain 
1 • u. The obvious left action of T on S* commutes with d. 

For any k > 0, let (S k )r be the module of T-coinvariants. This is the quotient 
of S k by the relations of the form 7 • u — u, where 7 G T, u G S k . This is also a 
complex with the induced boundary, which we denote by dr- It is known (cf. ]T5| ) 
that H u ~ k (T; C) is naturally isomorphic to H]<.((S*)t ® C). 

Let u = [t>i, . . . ,t> n +fc] be a fc-sharbly. Let ||u|| be 

Max I det(v h , . . .,v in )\, 



2 Since v G M(a) <^4> — v G M(cr), it is customary to write down only one member of the pair ±w. 
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where the maximum is taken over all n-fold subsets {ii, . . . , i n } C {1, . . . , n+k}. Note 
that this quantity is well-defined mod the relations in Definition |2.6| . We extend this 
notion to sharbly chains £ = J^n(u)u by setting ||£|| to be the maximum of ||u||, as 
u ranges over all sharblies in the support of £. We say that £ is reduced if ||£|| = 1. 
It is known (cf. PTfl ) that for T C SL 4 (Z), the group H 5 (Y; C) is spanned by reduced 
1-sharbly cycles. 

2.7. Since the generators of the sharbly complex are indexed by sets of primitive 
vectors, it is clear that there is a close relationship between S* and the chain complex 
associated to W, although of course S* is much bigger. Both complexes compute 
H*(Y;C). We refer to for a discussion of this, phrased in terms of the fan 
V . The main advantage of (S*)r is that it admits a Hecke action. Specifically, let 
£ = ^n(u)u be a sharbly cycle mod T, and consider the Hecke operator T(l,k) 
associated to the double coset YD(l,k)Y (cf. §|1.1|). Write 

YD(l,k)Y = ]jYg, 

g en 

a finite (disjoint) union. Then 

(2) *)(£)= Muk-u. 

Since Q (jL SL n (Z) in general, the Hecke- image of a reduced sharbly isn't usually 
reduced. 

3. Implementation details 

3.1. We state our results for general n as much as possible, though our main case of 
interest is n = 4. We have working programs for n < 4. Though we focus on SL ra (Z), 
analogous results hold for GL n (Z), and we have working programs for both SL and 
GL. 

Section |3] is very technical. The reader may wish to skip to §[| or §|5|. 

3.2. Let r = SL n (Z). Recall that Wt is the T-orbit of cells of type T in W. Let 

be a fixed representative cell in Wt- The stabilizer in f of ot is denoted f CTT , or Tt 
for short. This is a finite group that is straightforward to compute, since the minimal 
vectors M(ctt) are known. Our program maintains a database of the Tt- 
Standard facts about stabilizers give the following: 

Proposition 3.3. There is a one-to-one correspondence between cells a G Wt and 
cosets Y/Yt, given by ^ot jYt for any 7 G f such that a = ^ot- 

Under the smaller group Y, the f-orbit Wt breaks up into suborbits. If Y were a 
torsion-free group, Y\W would be a finite cell complex, its cells would be given exactly 
by the T-suborbits, and we could compute H l (Y\W) by the standard methods for cell 
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complexes. In our case, T is not torsion-free, but T\W can be thought of as a finite 
"orbifold cell" complex, whose elements are the quotients of cells by finite groups; 
the T-suborbits are in one-to-one correspondence with the orbifold cells. 

The goal of this subsection is to understand the T-suborbits in terms of the actions 
of the iV on finite projective spaces. By P n_1 = P n_1 (Z/./VZ), we mean the set of vec- 
tors (xi, . . . , x n ) G (Z/iVZ) 71 that are primitive in the sense that the ideal (x±, . . . , x n ) 
in Z/iVZ is (1), modulo the equivalence relation given by scalar multiplication by the 
units (Z/iVZ)* of Z/iVZ. When N is a prime, this is the usual projective space over 
the field of N elements. As usual, the equivalence class of the vector (x±, . . . ,x n ) is 
denoted a = [xi : • • • : x n }. We view these n-tuples as rows rather than columns; f 
acts on the right on P n_1 (Z/iVZ) in the obvious way. 

We define the bottom row map b : T — > P n_1 (Z/./VZ) as follows. For a matrix 7 G T, 
the bottom row of 7 is a primitive vector in Z n . Let 6(7) be the equivalence class of 
this image in P n ~ 1 (Z/A^Z). 

Lemma 3.4. The bottom row map b : T —>■ P™ _1 (Z/A^Z) induces a bijection between 
r\L and F n - l (Z/NZ), given by 

L7 t— > 6(7). 

The map is equivariant for the right action ofT. 

Proof. It is a standard fact that a vector in Z n is primitive if and only if it is the 
bottom row of some element of T. This implies b is surjective and that b _1 ([0 : ■ ■ • : 
: : 1]) = T. The rest is clear. □ 

We can now describe the T-orbits of cells in each Wj<. 

Proposition 3.5. The V -orbits of cells in Wt are in one-to-one correspondence with 
the orbits O of the right Vx-action on P n_1 . 

Proof. We have T\Wt = r\r/F T by Proposition K73L and this equals P n-1 /fY by 



Lemma □ 



The first step of our computer program is to determine, for each type T, the 
decomposition of P n_1 into right IY-orbits. Since we are primarily studying H^(W; C) 
for i = 5,6, we only need to work with the T representing cells of dimensions 4, 5 
and 6. 

We note the following fact, whose proof is immediate. 

Lemma 3.6. Let 70 G f , and let a = 6(70) in P n_1 . Then the stabilizer of a under 
the right action off is j^T'jo. 

Let a G W be a cell of type T, with a = ^o^t- Its stabilizer in T, denoted T a , is 
clearly 

(3) r CT = ( 7o f r7o - 1 )nr. 

Lemma |3.6| implies 
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Lemma 3.7. The group 7 1 r <T 7o is the subgroup of Tt that preserves a = 6(70) 
under the right action on P n_1 . 

3.8. In this subsection, we fix orientations on the cells a G W. It's necessary to 
be extremely careful — mistakes in orientation are easy to make and will ruin the 
computations. The price to pay is to sort through the details of the action of f . 

Recall that Tt is the stabilizer of o~t in T. For each T, there is an orientation 
character fV — > {±1} indicating whether or not 7 G Tt preserves the orientation on 
o~t- Our program stores the values of these characters along with Tt- Let be the 
subgroup of Tt where the orientation is +1. 

Remark 3.9. To compute the value of the character at 7, we determine (i) how 7 acts 
on the orientation of the cone C, and divide by (ii) how 7 acts on the orientation of 
the Voronoi cone dual to cy. Dividing works because C is locally the direct product 
of the cell and its dual Voronoi cone. As for (i), every element of T acts by +1 on 
the orientation of C, since T is a subgroup of the connected group SL n (R). In this 
paper, where n = 4 and dim a > 0, it turns out that all the dual Voronoi cones 
are simplicial; the sign in (ii) is the sign of the permutation that 7 effects on the 
bounding rays q(v) of the cone, which is easily computed. 

Let O be a right IY-orbit in P™" 1 . We call O non-orientable if for some (which 
implies every) a G O, there exists 7 G Tt n LJ with a7 = a. Otherwise, we call O 
orientable. These notions depend on T, though we usually leave T out of the notation. 

If O is orientable, fix some ao G O. Define the orientation number of a G O to 
be +1 (resp., —1) according as a = a 7 for some 7 G f J (resp., 7 G IV N fj). 
The orientation number is well-defined precisely because O is orientable. Again, the 
notions depend on the choice of a , though we leave a out of the notation. 

Let 7 G IV, and let p be any cell of W with any given orientation. Since 7 acts by 
diffeomorphisms on C, it carries the orientation on p to some orientation on the cell 
7p. We write 

(4) (7)*(P) 

to denote jp together with this orientation. Clearly (7)* is functorial, and preserves 
the relative orientation of p, t whenever r is a co dimension-one face of p. 

Once and for all, fix orientations on the standard cells ot- We can now put 
orientations on all the cells of W. 

Definition 3.10. Let a be a cell in Wt with a = 7o0t. Let O be the right T^-orbit 
in P n_1 corresponding to a as in Proposition |3.5| . If O is orientable, we give a the 
orientation 

(5) (orientation number of a) ■ (70)* (or). 
If O is non-orientable, we give a an arbitrary orientation. 
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Proposition 3.11. The quantity in (|5|) is well-defined. 

Proof. Let a, 70, O, and a be as in Definition |3.10| , with O assumed orientable. 
Assume a = 71 07 as well as 7007. Let a x = 6(71). Then 7f 1 7o G Fy. By definition 
of f J, (7o)*(o"t) = (7i)*(cr r ) if and only if 7-f 1 7 G f J. On the other hand, ai7f/ 70 = 
[0 : • ■ • : : l]7o = a, so a and ai have the same orientation number if and only if 
7r 1 7oGf+ ' □ 

We must understand how T and T act on the orientations we have just chosen. 



The following lemma is immediate from (f|) and Definition p. 10 



Lemma 3.12. If a = 70 cr^ corresponds to an orientable O, thenjo carries the chosen 
orientation of o~t to the chosen orientation of a times the orientation number of a. 

Here is a more general statement. 

Proposition 3.13. Let a = 7007- and 0\ = jia, for some 70,71 G T. Let a = 
b(7o) and &i = 6(7170). Let 0\,0 be the right T T -orbits in P"- 1 containing a!,a ; 
assume these orbits are both orientable. Then 71 carries a to U\ while multiplying 
the orientations by 

(6) (orien. number of a) ■ (orien. number o/ai). 

Proof. Apply Lemma |3.12| twice. □ 



Fortunately, @ becomes trivial when we consider T as opposed to V. 

Proposition 3.14. Let o\ = , -fio~ for 7! £ T. Let a 1; a, Oi,0 be as in Proposi- 
tion \3.l% both orbits being assumed orientable. Then carries a to o\ with orien- 



tations matching. 

Proof. We have [0 : • • • : : 1)71 = [0 : • • • : : 1] by the definition of T. Hence 
ai = ([0 : • • • : : 1] 71)70 = [0 : ■ ■ • : : l]7o = a. Thus the expression in (||) is a 
square either of +1 or of —1. □ 

3.15. We compute the equivariant cohomology H£(W; C) using a spectral sequence, 
following the exposition in [13|, VII. 7-8]. (The spectral sequence there is for equivari- 
ant homology; we make the appropriate modifications for cohomology.) 

Let a be any cell in the well-rounded retract W. Recall that r CT is the stabilizer 
of er in T. Let be the T^-module where 7 G T a acts by +1 if 7 preserves the 
orientation of a and by —1 if it does not. 

For each i, let Wu\ be a fixed set of representatives of the T-orbits of the cells in 
W of dimension i. The E% term of the spectral sequence is 

(7) E^= & p.; C). 



(0 
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These cells o (omicron) are in one-to-one correspondence with the O of Proposi- 



tion |3.5| , as T runs through the types of cells of dimension i. Because C is a field of 
characteristic zero, all the terms in (|7]) vanish when j ^ 0. In particular, the spectral 
sequence collapses at E 2 . The term H (T o ; C D ) is the subset of r o -invariants in the 
module C Q . 

Proposition 3.16. For any cell a, let T be the type of a, and let O be the right 
Tx-orbit in P n_1 corresponding to a as in Proposition |<?. & Then H°(T a ; C) is C if O 
is orientable, and is if O is non-orientable. 

Proof. When O is orientable, this follows from Proposition [3.14| , merely because T a C 



T. Now assume O is non-orientable. Let a = 70 or, with a = b(7o). As we have said 
above, there is some 71 G T T \ PJ with sr/i = a. By (2.2), T a = (7ofT7o ) R T. The 
element ■JoJi^q 1 is i n r by Lemma 2, so it is in r CT ; clearly it carries a to itself while 
reversing orientation. Hence r CT acts non-trivially on C a , meaning H°(T a ; C CT ) = 
0. □ 

Remark 3.17. The proposition shows that the T-orbits of cells coming from non- 
orientable O contribute nothing to our spectral sequence. We ignore these objects 
for the rest of the computation, tacitly assuming that all O's mentioned from now 
on are orientable. 

To summarize: 

Proposition 3.18. The E\'° term of the equivariant cohomology spectral sequence 
for ifp(W;C) is a direct sum ® C ; where o runs through a set of i- cells in one-to- 
one correspondence with the orientable right Tx-orbits in P n_1 , for all types T of cells 
of dimension i. 

3.19. In § |3.21| , we will describe the boundary maps di of the spectral sequence. 



These are the only differentials we need consider, since the sequence collapses at E 2 . 
In this subsection, we give some details concerning how the cells meet at their bound- 
aries. 

As usual, a facet of a cell a is any face of a of codimension one. Let be the set 
of facets of a in W. 

We will need to understand how J a breaks up into orbits under the action of T a 
(and to choose a set J' a of representatives of these orbits in 3 r a ). We do this in 
Proposition |3.2U| below. We will start by determining 5F CTT for the standard cells o~t 



in a form suited to our computation. We will then determine 3^ for any a. 

We make two conventions, (i) If types T and T 1 occur in the same discussion, it 
is assumed that a cell of type T has at least some cells of type T' as facets, (ii) If 
a = xo~t for some x G F, we identify a with the coset xTt as in Proposition |3.3| . In 
the expressions of the form 

U(*)> 
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the (*) will be a finite union of cosets, say X\Tt' ]J • • • ]J x^Yt 1 for Xi,...,Xk £ 
r. It is understood that (*) corresponds to the set of cells Xi(Jt', ■ ■ ■ ,Xk<JT', as in 
Proposition |Q[ Even when (*) is a more complicated object, like a double coset, its 
meaning is that one should decompose it into single cosets by choosing appropriate 
representatives (which will be the xi). 

The boundary of <jt is a union of (the closures of) various cells of type T'\ 



= \Jf3 T >,iT TI \\ ■ ■ ■ Y[p T , k T T , 



for some ^T',i £ r. However, (K|) is invariant under the left action of cr^'s stabilizer 
Ft- Hence there must be finitely many a(T,T',t) G T such that 



(9) 



3v T = [JrYa(T,T',i)rV Y\_ ' ' ' TT rrQ(r,r',fc)iV- 



We have computed the oc(t t',i) by hand for n < 4. In our cases of interest (n = 4, 
dima > 4), one finds there is only one l — that is, the right-hand side of (|D|) is actually 
just one double coset. In fact, we find that we may take ce(T,T',t) to be the identity 
except in one case, where (T, T") = (5b, 4b) in the notation of PT| and 



«5b,4b,l 



/ 1 \ 

1 

-10 10 

\ -1 1 ) 



From now on, we will write a for ctrrT',i)i the dependence on T and T' being under- 
stood. 

To find an expression for J a for a = 7oC"t, we multiply (pi) by 7 to obtain 



(10) 



= [J 7orTaIV 



To exhibit this as a union of cosets of the form XjIV, we rewrite it as 

(11) 9 a = \Jj a(a- 1 f T a)f T , 

(12) = I^J r y a(a~ 1 T T a/ (a _1 f T a fl f r'))f t>, 



The last formula (for each T") is a disjoint union of single cosets, in one-to-one 
correspondence with a set of representatives of the cosets in a~ l Y T a/ (a~ 1 r r a;rir T /). 
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It is a matter of formal manipulation to get an expression for Y a \J a : 

(13) rv\^ = |J(7of no 1 n r)\ 7o «(a- 1 f T a/(a-T T a n r T ,))T T , 

(14) = |J 7o (r T n 7o ~ 1 r7o)\«(«- 1 r T a/(a- 1 r T a n r T ,))r T , 

T' 

(15) =1 j7o«[a- 1 (r T n7 ( r 1 r7 )a\a- 1 f T a/(a- 1 f T anf r /)]f T / 

1 A C 



where for each T' the right-hand side of fllS]) is again expressed as a disjoint union of 
fV-cosets, in one-to-one correspondence with a set of representatives of the double 
coset expression in square brackets. 

We must interpret ([TT]) (|T2D in terms of the fY'-orbits in P n_1 . As usual, let 
O be the right IY-orbit corresponding to a, with a = b( 7 o). Let c = b(joa), so 



that c = aa. Now the quantity A in (15) is exactly the subgroup of a 1 Ytcv that 



preserves c. Hence any set of coset representatives for B in (|T5|) is a set of matrices 
whose bottom rows, in P n_1 , are the members of the right (a _1 f ^a)-orbit of c. 
Equivalently, any set of representatives for B is a set of matrices whose bottom rows 
are exactly the members of O ■ a. The group denoted C in (|T5| ) acts on this orbit 
O ■ a, decomposing it into suborbits; for given T', the disjoint f ^/-cosets in (pTT|) — (pT2|) 
are in one-to-one correspondence with these suborbits. We summarize this result as 
follows: 



Proposition 3.20. Let T , T' and a be as introduced in this subsection. Let o be 
a cell of type T, represented by the T T -orbit O in the manner of § |ff. | . Decompose 
the orbit O ■ a into its suborbits 0\, ... ,Ok under the group C = a^fVa D Fy. Let 
a.j G Oj. Let 7 j G f be chosen so that b( 7j -) = a,. (We may, in fact, take 7j - to be of 
the form 707^0; for some 7,- G Yt-) Then the union over all T' of the cells 

{ 7 10T', • • • ,7fc°"T'} 

is a set J' a of representatives for Y a \S a . 

For all pairs T,T', we store the intersection C = a~ l Yxct n Yt' in our program. 



3.21. We can now determine the boundary maps d\ in the spectral sequence. Recall 
that W(i) is a fixed set of representatives of the T-orbits of the cells in W of dimen- 
sion i. Superseding the use of o in § [3.151 , we let a run through Wfy+i), and let r run 
through the set 3^ of representatives of the facets of a. To compute H^(W;C) for 
subgroups of SL 4 (Z) for i = 5,6, we must work out the map d % {° for i = 4 and 5. 
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We follow ||T3|, VII. 8], taking the dual to turn homology into cohomology. The map 

(16) df°: H°(T T0 ;C T0 ) - H°(r a ;C a ) 

T £W {i) rt%l) 

is a sum of terms, one for each pair (a, r); here r is the fixed representative in Wu\ 
that is T-equivalent to r. From now on, we focus on a single pair a, r. Call their 
types T and T', respectively. Let O be the right f ^-orbit in P™" 1 associated to a by 
Proposition [375]. If a = 70 or, let a = K7o) e Let To-r = To- H T r . 
As in pT3] , p. 176], the term for (a, r) in g^' is the composition 

(17) H°(T T0 ; C T0 ) tf°(r r ; C T ) tf°(r CTT ; C CT ) ^°(r ff ; C CT ) . 

We now give the definition of these maps and show how to compute them. Note 
that all the coefficient modules are copies of C on which the groups act trivially. 

The map t aT is the transfer map C — > C given by multiplication by the scalar 
[IV : T aT \. We evaluate this scalar as follows. First of all, T a = ^oFtTcT 1 ) H T, which 
has the same cardinality as T T D (•yQ 1 Tj ). The latter is the subgroup of T T that 
fixes a. Thus 

(is) #(r^ #(Ft) #(Ft) 



#(r T -orbit of a) #(0) ' 

The numerator is known because we stored Ft- The denominator is easily recovered 
from the computer's lists of orbits. 

We now evaluate the order of F ar . As before, F a = 7 f tIo 1 nf. By Proposi- 
tion |3.20| , t = 7o7ao"T' for some 7 G Ft- Hence 

T r = (( 7 o7a)f T ,(a- 1 7~ 1 7o _1 )) n T. 

Writing out F aT = Y a r\Y T and conjugating by (707a), we find that has the same 
cardinality as 

(19) (a~ l F T a n t T >) n(a _1 7 _1 7o" 1 r7o7a) 

V V ' 

D 

The group D is the same as the group C from (|T5|), and is one of the groups we've 
stored. And the whole group in fll9"|) is exactly the subgroup of D that fixes the point 
ay of Proposition |3.20| . So 

(20) Mr x #("~ 1 rranr g v) 
(20) #(r ffT ) _ _ . 

Again, the numerator is known from what's stored, and the denominator is easy to 
evaluate. 
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We have 

(21) taT = m' 

As a consistency check, the program signals an error if the computed value of ([H]) is 
not an integer. 

Brown's u aT : C — > C is the composition 

H°(T T ; C T ) -> i?°(r CTT ; C T ) -> tf°(r aT ; <C CT ). 

The first arrow is induced from the inclusion T aT ■=— > T T , and is easily seen to be 
the identity. The second arrow is induced by the r (TT -map d aT : C a — ► C T , namely 
the (cr, r)-component of the cellular boundary operator on W. Let [a : r] be ±1 



depending on whether the orientation on a from Definition 3.10 does or does not 



induce the orientation on the facet r from Definition 3.10 . Then d aT is the map 
C — ► C given by the scalar [cr : r]. Thus m (7T = [a : r]. 

To evaluate [a : r] , we reduce the problem to evaluating the boundary operator on 
a small list of pairs of cells. Write r = •y Q ^aaT'- Then a = 707 or, since 7 G IV. If 
x G P n_1 is part of any oriented IV-orbit (for any T), write sgn T (x) for the orientation 
number of x with respect to IV. By Definition |3.10| , a and r — together with their 
standard orientations — are given as follows: 

(22) a = (707) *M ■ sgn T (b(7o7)) 

(23) t = (707a) *( cr T / ) sgn T , ( b (707a) ) • 
Because (...)* is functorial, 

t = (7o7)*(a)*((7 T /) • sgn T ,(b(7o7«)). 
Since (...)* preserves the [ : ] relation, we may cancel out (707) *'s, obtaining 

(24) [<t:t] = sgn T ( b (707) ) ■ sgn T , ( b (707a) ) ■ [a T ■ (a) # (a T r)]. 

E F 

We evaluate each factor in (fM|) in turn. In F, b(7o7a) is the point of Propo- 
sition |3.20| , and sgn T ,(aj) is simply the orientation number of in its r^z-orbit. 
Similarly, E is the orientation number of 6(707), which is a T^-translate of a, in its 
IV-orbit. In practice, though, we do not know 7 or 6(707) explicitly, so it is easier 
to evaluate E by another method. Note that b(7o7) ■ a = 6(7070) = a^, an element 
of the (a _1 f r a)-orbit O ■ a of Proposition |3.20| . It is easy to compute all the points 
in O ■ a, using the group a _1 IVa (which was stored). For some x, we will have the 
equation a^ = x • a. Then E will be the orientation number for this x. 

The quantities \<tt '■ ( Q; )*( cr r')] are evaluated by hand, for all T, T' occurring in (0). 
Several issues arise. First, we must find [<jt '■ &t'] whenever the standard cell <jt< is 
a facet of cy. No matter how we choose the orientations at the start, it is in general 
impossible to arrange our choices so that all the relative orientations are positive. As 
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a general illustration, if a 0-cell, two 1-cells, and a 2-cell meet locally in a picture 
like the first quadrant of M 2 , and if both 1-cells are oriented to point away from the 
origin, then the pair (2-cell, x-axis) must have relative orientation opposite to that 
of the pair (2-cell, y-axis), no matter how we orient the 2-cell. In the SL 4 (Z) case, 
one can draw a schematic picture of how the standard cells meet and can read off all 
the relative orientations [<jt : &t']- 

Second, we must compute (a)*(<TT') when a isn't the identity, knowing the orien- 
tation on <jt< • This involves the same techniques as in Remark [3.9| . 

Third, we must find some facet v of ot whose orientation is known, and must 
compare [cr T : v] to [a T : (a)*(t7y/)]. In practice, if dim cr-r = k, this means finding 
a chain of (k — l)-cells between v and (ol)*{o~t') such that consecutive members of 
the chain meet in faces of dimension k — 2, and comparing the orientations of v and 
(a)*(<7y/) across the (k — 2)-faces. 

Finally, we must compute v T , but this is easy. It is induced by the conjugation 



action of the element of T that carries r to tq. But by Proposition |3.14| , any element 
of T preserves the orientations on the cells. Hence one finds v T = +1. 

4. Hecke operators 

4.1. We identify the cochain complex in flTop with the complex of cellular cochains 
on W by identifying an i-cell a with a generator of H°(T a ] C CT ), taking care to make 
the signs match. Formulas like ^n(cr)cr will denote the corresponding cocycles in 
either complex, and will be referred to as W -cocycles 

Let (3 G H 5 (T;C) be a class, and let u = ^tz(ct)<7 be a representative for (3 in 
terms of the previous paragraph. Let T(l, k) be a Hecke operator. To compute the 
action of T(l, k) on u, we do the following: 

1. Convert u to a reduced 1-sharbly cycle £, and then compute T(l, k)(£) using (^j) 
in 

2. Use the algorithm from [IS| to write T(l,k)(£) as a sum of reduced 1-sharbly 



cycles. 

3. Convert these reduced 1-sharbly cycles to PU-cocycles. 

Step is described in detail in |[L8|| , and we refer the reader to that article. In this 
section, we focus on steps [I] and |] in the context of §|3|. 
We begin with a definition from WF 



Definition 4.2. [T^ Definition 5.3] Let u be a basis element of the fc-sharblies Sk- 
Then a lift for u is an n x (n + k) integral matrix M with primitive columns such 
that [Mi, . . . , M n+ k] = u, where Mj is the ith column of M. 

Modulo the action of GL 4 (Z), there is only one orbit of reduced basis O-sharblies 
and only three orbits of reduced basis 1-sharblies. The identity matrix serves as a 
lift for a member of the first orbit, and lifts representing elements of the latter three 
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orbits are 

/ 1 





1 \ 



/ 



/ 1 






i \ 
i 
i 



and 



/ 1 





i \ 
i 





We call these the standard 0- and 1-sharblies. The sets of primitive vectors indexing 
the standard 6 and 5-cells in W coincide with the sets of column vectors of these 
matrices. By abuse of language we will speak of the "standard sharbly of type T," 
and will use the notation u^. 

4.3. Given a sharbly cycle £, we denote by supp£ the support of £. Suppose that 
T is torsion-free. Then according to [P|, a 1-sharbly cycle £ mod V with coefficients 
in a ring R can be encoded by a collection of 4-tuples (u, n(u), {v}, {L(v)}), where 

1. u G supp £, 

2. n(u) G R, 

3. {v} = supp<9u, and 

4. {L(v)} is a T-equivariant set of lifts for {v}. 

The T-equivariance condition in |4] is the following. Suppose that for u, u' G supp £ 
there exist v G supp(<9u) and v' G supp(<9u') such that v = 7 • v' for some 7 G T. 
Then we require L(y) = r yL{y r ). 

In the case under study, Y is not torsion-free, and the above data needs to be 
modified. Suppose that a 0-sharbly v G supp du has a nontrivial stabilizer T(v) C T, 
and let m be any lift of v. Then in the cycle £ we replace n(u)u with 



n u 



^ #r(v) 

7 er(v) ^ v y 



u 



1 ! 



where u 7 has the same data as u, except that we give v the lift 7m. (Note that this 
is possible in our case since the coefficient ring R = C is divisible.) 



4.4. Now we describe how to construct the data in §fO]to produce a 1-sharbly chain 
£ corresponding to the W-cocyde u. There are two steps. 

First, choose a such that n(a) 7^ in u. According to Proposition |373|, the 5- 
cell a G Y\W is encoded as a coset [7IY], where 7 G V and Y T is the stabilizer of 
the standard cell of type T. Moreover, the coset [7IY] is encoded as the orbit O, 
which can be regarded a set of triples {(a, ±l,7 a )}, where a G P 3 (Z/iVZ), 7 a G T T , 
and ±1 is the orientation number (§ |3.8|) . From O we arbitrarily choose a triple 
with orientation number 1, and then using Hermite normal form construct a matrix 
7 G SIj4(Z) with bottom row equal to a. Then the contribution of o to the 1-sharbly 
chain is 



n(cr)7U T . 
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We do this for all of supp-u and sum to produce £. For each u G supp£, we write 
11(7) if we want to indicate the element 7 used in the construction of u. 

4.5. At this stage, we have a 1-sharbly chain, and we need to choose lifts to reflect 
the cycle structure of £. This we do as follows. In the spirit of § 3.19| , for each type T 



we choose a set of matrices VIt C SL^fZ) such that d: Si — > So can be written as 



d: u T 1 — ► 



wv, 

where v is the standard basis 0-sharbly. Note the absence of signs in this map — the 
signs in the boundary map in Definition |2.6| have been encoded in the lu's, which may 
nontrivially permute the column vectors of v. 
Form the 0-sharbly chain 

(25) Yl W 

u(7),aien T 

where v is the standard 0-sharbly, we sum over all 11(7) G supp £, and T is the type 
of 11(7). (In ( pB"l) we have abbreviated a(T,T',i) to aii, since T is determined by u(7), 
and there is only one type of reduced 0-sharbly mod SL 4 (Z).) Note that this sum is 
in if?*, not [S*)r] the only relations we apply are those in the sharbly complex. 

4.6. After summing, we find that some O-sharblies cancel, and some remain. For 
those that canceled, we can choose any lifts we like, as long as we choose the same 
lifts for all terms that cancel each other. 

The remaining O-sharblies form a chain rj that vanishes in (So)r <8> C, and we must 
choose nontrivial lifts for them. To do this, first arbitrarily choose lifts for each 0- 



sharbly in supp rj. The data we computed in § |3.2| allows us to easily compute the 
distinct orbits of Y in supp?]. We do this and order each orbit. 

Suppose v is the first 0-sharbly in one of these orbits, and that it corresponds to 
the triple (ao, ±l,7 ao ). Let vo(7) be its lift. Then if v is any other 0-sharbly in vo's 
orbit, corresponding to the triple (a, ±l,7 a ), we replace the lift v(7) by 

V(7) < V (7)7a Sao 

After these lifts are constructed, the cycle £ is ready for input in the Hecke operator 
program. 

4.7. Upon completion, the Hecke operator program returns a reduced 1-sharbly 
cycle, which we must convert to a PU-cocycle. So let £ = ^n(u)u be a reduced 
1-sharbly cycle, and let u G supp£. First we determine which of the three types T 
of standard reduced 1-sharblies u has. Then we must find a matrix 7(11) G SL 4 (Z) 
such that 

7(11) • u T = u, 
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where u-p is the standard 1-sharbly of type T. This is straightforward, although one 
must be careful to incorporate the orientation number of 7(11). 

In practice, the main step is the following. Let u be a reduced 1-sharbly basis 
element with lift M. We choose a nonsingular 4x4 minor m of M and construct 
m . Then m _1 • u will be a 1-sharbly with lift m _1 M, and will be standard except 
possibly for one column vector. By multiplying m~ l M on the left by elements of the 
stabilizer of the standard 0-sharbly, we can eventually produce the standard 1-sharbly 
with same type as u. This allows us to construct 7(11). 

5. Numerical results 

5.1. In this section we present numerical data from our experiments. As mentioned 
in §[[], to avoid floating point problems with C-coefficients we usually work with 
F = F31991, the finite field with 31991 elements, and in some cases work with Z or Q. 
The computations were carried out on a variety of Unix machines at Columbia and 
Oklahoma State. The code for the cohomology of W (§13) was written in Common 
Lisp. The Hecke operator code (§|j) was written in C++, and used the LiDIA library 
||14||. Perl scripts patched together the outputs of the various programs and produced 



the tables in 55.5 



5.2. We first describe how we performed linear algebra on the large sparse matrices 



that arise in our computations. Fix the level N. We use the notation of § |3.21 
working over C at first, and letting Vu) be the domain of the map d l f in ([!§). To 
compute if 5 (T;C), we must find the kernel of d^° modulo the image of 4'°. We 
prefer to find the kernel of a single matrix. Regard the d\ as matrices acting on the 
left on column vectors, and let D% be the matrix where g^' is stacked on top of the 
transpose of d{ : 



IDs 



d{ 



This matrix defines a map V( 5 ) — > V( 6 ) © V( 4 ), where we use the standard inner product 



to identify the transpose of d±'° with its adjoint. The kernel of this map is the space 
of harmonic 5-cocycles; it is isomorphic to if 5 (T;C). Let s = dimV^, the number 
of columns of ®z- 

The matrix Dz has coefficients in Z. For any ring R, set Qr = £>z <8> R, and set 
D = 2)f- In the tables in § |5.4| , the value of "rank" we report in the rows labeled R 
is a number almost certainly equal to dimkerDj?, whose computation is explained 
below. This number is also almost certainly equal to dim H 5 (T; C). 

5.3. When 25 is very large, we could not have found its kernel without a sparse ver- 
sion of the Lanczos algorithm. This algorithm is usually used with real or complex 
matrices, particularly for eigenvalue problems. Following ideas in we translated 



it into the mod-p setting. Let (£ = D ■ D, a symmetrized version of D. We choose 
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a random non-zero seed vector v with coefficients in F and consider the sequence 
v, (Ev, (E 2 v, .... The Lanczos algorithm shows us how to compute not this sequence, 
but the sequence v = q , qi, ... resulting from it by the Gram-Schmidt orthogo- 
nalization process. We perform the Gram-Schmidt process mod p in the naive way, 
using Yl x jVj f° r the inner product. This means we don't have the usual guarantee 
that ^2 x "j wm be non-zero when (xi, X2, ■ ■ ■ ) 7^ (0, 0, . . . ). If the inner product is 
ever for non-trivial (xj), we simply abort and choose another random seed v; even 
for large D, these aborts happen less than half the time. 

The strength of the algorithm is that the RAM only has to hold the sparse matrix 
D and a few vectors of storage. It does not have to hold (£, which is dense in general. 
The q's form a dense matrix, but they may be stored on the disk, not in RAM. (In 
our implementation, D tr was stored in RAM along with £).) 

Let k be the largest value for which the set {q , . . . , q s ~k} is linearly dependent. 
Reading the q's back in from disk, the algorithm allows us to backsolve for a non- 
zero vector y G ker (£. What we want is an element of ker D. A priori, we only know 
ker (£ D ker D, and we will see that the containment is not always an equality (though 
it would be over 1 or C). However, by checking 2) ■ y = directly, we always find in 
practice that the y we compute lie in ker 2). 

Thus each successful run of our algorithm produces one kernel vector for D. It 
also produces k. One can easily show k > dim ker (E. We run Lanczos up to 30 or 
40 times with different random seeds, and we find k is independent of the random 
seed used.0 We conclude that dim ker € = k; though we have not proved this, the 
computational evidence seems conclusive. It is clear that dim ker (£ > dim ker D. 

To find a basis of ker 2D, we run Lanczos many times until we have a set S of k + 10 
elements of ker D. We use mod p Gram-Schmidt on subsets S' C S to find maximal 
linearly independent subsets of S. We start with several different <S"'s. In all the 
cases we checked, we found that maximal linearly indepdent sets in S had a common 
cardinality k', that any fc'-element subset of S was linearly independent, and that 
any subset of S with more than kl elements was dependent. We conclude that k' is 
the value of dim ker 2); again, the computational evidence is convincing, though all 
we have proved is k' < dim ker D. In the rows marked "Lanczos" in Table [I], what 
we report as "rank" is k! . In the rows marked Z or Q, what we report (namely k') is 
provably the rank of H 5 (T, C). We do find k! < k sometimes in practice. 

5.4. In Table p], we give the results of our Betti number computations. Gaufi means 
that ordinary Gaussian elimination was used to find the kernel of D, and Lanczos 



3 For one level, random seeds produced a certain value of k, while one random seed produced a 
value that was greater by 1. The Lanczos method behaves this way when E = {(5v, € 2 v, (£ 3 v, . . . } 
does not span the image of £, but only a proper subspace of the image. Since v is chosen randomly, 
it is extremely rare for E to span less than the full image; our data bears this statement out. 
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means the algorithm of § [5.3| was used. The entries marked with * are those for which 
k — k'. 



I ,fMrr>l 




ran \c 


1 Cllldl JA.O 






ran [r 


t'PTI 1 Q tItC 
I CLllClI rvo 


1 1 
11 


L 


r> 

Z 


(oraUD 


on 

33 


TI? 

It 1 


10 


Lanczos 


1 o 

12 


TT7' 

If 


u 




34 


If 


12* 




13 


iL 


l 




35 


TT7' 

If 


7* 


57 


14 


Jr 1 


2 




3d 


If 


24 


57 


15 


TI? 


Z 




37 


TT? 
If 


8 


GauD 


16 


TI? 
If 


3 




38 


TT? 
If 


1 A ^ 

14 


Lanczos 


17 


IL 


3 




on 

39 


TI? 
If 


1 n^ 

10 




lo 


TI? 
If 


9 




A C\ 

40 


TI? 
If 


9 




19 


IL 


3 




A 1 

41 


TI? 
If 


n^ 

9 




20 


TT? 
If 


r> 

Z 




A O 

4z 


TI? 
If 


17 




zl 


TT? 
If 


3 




A O 

4o 


TI? 
If 


10 




99 
zz 


F 

Jr 


7 






F 

Jr 


lo 




23 


z 


5 




45 


F 


27 




24 


F 


2 




46 


F 


19* 




25 


F 


7 




47 


F 


11* 




26 


F 


7 




48 


F 


26 




27 


F 


12 




49 


F 


33* 




28 


F 


7 




50 


F 


34 




29 


F,Q 


6 




51 


F 


19* 




30 


F 


8 


Lanczos 


52 


F 


21* 




31 


F 


6 


Gau8 


53 


F 


17* 




32 


F 


12 


Lanczos 











TABLE 1. Probable Betti numbers for H 5 (T (N); C). 



5.5. Next, we present the Hecke data we computed. These computations are much 
more arduous than computing Betti numbers, and grow in complexity very fast as 
a function of the number of cells of the retract W mod T. Hence we were able to 
compute only a few Hecke operators, usually only T(2, *) and T(3, *). Beyond level 
20, it becomes infeasible to compute T(3, *); hence most of our data at large levels 
is only for T(2, *). Happily this is usually sufficient to guess persuasively what is 
happening with the cohomology. 

We use the following conventions and abbreviations. All polynomials should be 
considered as elements of F[X], where the bar denotes algebraic closure. We denote 
the p-adic cyclotomic character of Gq by e, so that e(Frobz) = I for any I ^ p. 
The symbol Ha (resp. lib, IV) denotes a Galois representation of the form eVfc © 



COHOMOLOGY OF CONGRUENCE SUBGROUPS OF SL 4 (Z) 



21 



e b © e c , where (k, a, 6, c) is (2, 0, 2, 3) (resp. (2, 2, 0, 1), (4, 0, 1, 2)), if a k is the Galois 
representation associated to a weight k classical holomorphic cuspidal newform / G 
S^ ew (N'), where N' divides N. The same symbol prefixed with an "E" denotes a 
Galois representation of the same form except that a k is the Galois representation 
attached to an Eisenstein series of weight k for Ti(N'). 

The individual tables are organized as follows. For each level we give the rank, as 
defined at the end of § |5.3| . Then each block gives data for the Hecke eigenspaces. 
The first column gives the type Galois representation seemingly attached to this 
eigenspace. The second column gives the dimension of this eigenspace. The third 
column gives the index I of the Hecke operator, and the fourth column gives the 
corresponding factored Hecke polynomial. This is the polynomial defined on the 
right of ([]]) in Definition L2; it succinctly encodes the Hecke action on any vector 
in the eigenspace. At the bottom of each table, we indicate the Q-splitting of the 
spaces of newforms S& = Sl° w (N) under the action of the Hecke operators; this data 



is from [23]. After each table, we comment of the eigenclasses. 



Level 11. rank = 2. 



Ila 


1 


T 2 


(1 


- 4X)(1 - 8X)(1 + 2X + 2X 2 ) 






n 


(1 


- 9X)(1 - 27X)(1 + X + 3X 2 ) 






T 5 


(1 


- 25X)(1 - 125X)(1 - X + 5X 2 ) 






T 7 


(1 


- 49X)(1 - 343X)(1 + 2X + IX 2 ) 


lib 


1 


T 2 


(1 


- X){\ - 2A)(1 + 8X + 32X 2 ) 






n 


(1 


-X)(l- 3A)(1 + 9X + 243X 2 ) 






T 5 


(1 


- X)(l - 5A)(1 - 25A + 3125A 2 ) 






T 7 


(1 


-X)(l- 7X)(1 + 98X + 16807A 2 ) 



dimg 2 (ll) = 1, dimg 4 (H) = 2 



The weight 4 newform doesn't lift. 



Level 13. rank = 1. 


IV 


1 




(1 - 2X)(1 - 4X)(1 + 5X + 8A 2 ) 






n 


(1 - 3A)(1 - 9A)(1 + 7X + 27 X 2 ) 






T 5 


(1 - 5A)(1 - 25X)(1 + 7X + 125X 2 ) 






T 7 


(1 - 7X){1 - 49X)(1 + 13X + 343A 2 ) 


dimS 2 (13) = 0, dimS 4 (13) = 


= 1 + 2 



Only the rational weight 4 newform lifts. 



Level 14. rank 


= 2. 






Ila 


1 


n 


(1 - 9A)(1 


-27X)(1 + 2A + 3A 2 ) 


lib 


1 


T 3 


(1-A)(1- 


3A)(1 + 18A + 243X 2 ) 


dim S 2 (14) = 1, 


dim S 4 (14) = 


= 1 + 1 
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No weight four newforms lift. 



Level 15 


. rank = 2. 




Ila 


1 


T 2 


(1 - 4X)(1 -8X)(1 + X + 2X 2 ) 


lib 


1 


T 2 


(1 - X)(l - 2X)(1 + AX + 32X 2 ) 


dimS 2 (15) = 1, dimS 4 (15) = 


= 1 + 1 



No weight four newforms lift. 



Level 17. rank = 2. 


Ila 


1 




(1 - 4X)(1 -8X)(1 + X + 2X 2 ) 






n 


(1-9X)(1-27X)(1 + 3X 2 ) 






T 5 


(1 - 25X)(1 - 125X)(1 + 2X + 5X 2 ) 






T 7 


(1 - 49X)(1 - 343X)(1 - AX + IX 2 ) 


lib 


1 


T 2 


(1 -X){1- 2X){1 + 4X + 32X 2 ) 






n 


(1 -X)(1 - 3X)(1 + 243X 2 ) 






T 5 


(1 - X)(l - 5X)(1 + 50X + 3125X 2 ) 






T 7 


(1 -X)(l- 7X)(1 - 196X + 16807X 2 ) 


dimS 2 (17) = 1, dimS 4 (17) = 


= 1 + 3 



No weight four newforms lift. 



Level 19. rank = 3. 


Ila 


1 


T 2 
T 3 
T 5 


(1 -4X)(1-8X)(1 + 2X 2 ) 

(1 - 9X)(1 - 27X)(1 + 2X + 3X 2 ) 

(1 - 25X)(1 - 125X)(1 - 3X + 5X 2 ) 


lib 


1 


T 2 
T 3 
T 5 


(1-X)(1-2X)(1 + 32X 2 ) 

(1 - X)(l - 3X)(1 + 18X + 243X 2 ) 

(1 - X)(l - 5X)(1 - 75X + 3125X 2 ) 


IV 


1 


T 2 
T 3 
T 5 


(1 - 2X)(1 - 4X)(1 + 3X + 8X 2 ) 
(1 - 3X)(1 - 9X)(1 + 5X + 27X 2 ) 
(1 - 5X)(1 - 25X)(1 + 12X + 125X 2 ) 


dimS 2 (19) = 1, dimS 4 (19) = 1 + 3 



The rational weight four newform lifts. 



Level 20. rank = 2. 


Ila 


1 


n 


(1 


- 9X)(1 - 27X)(1 + 2X + 3X 2 ) 


lib 


1 




(1 


- X)(l - 3X)(1 + 18X + 243X 2 ) 


dimS 2 (20) = 1, dimS 4 (20) -- 


= 1 





The rational weight four newform doesn't lift. 
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Level 21. rank = 3. 


Ila 


1 


T 2 


(1 - 4X)(1 -8X)(1 + X + 2X 2 ) 


lib 


1 




(1 - X)(l - 2X)(1 + 4X + 32X 2 ) 


IV 


1 


T 2 


(1 - 2X)(l - 4X)(1 + 3X + 8X 2 ) 


dim ^(21) = 1, dimS 4 (21) = 1 + 1 + 2 



Of the weight four rational newforms, only one lifts. 



Level 23. rank = 5. 


Ila 


2 


T 2 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 


lib 


2 




(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 + 2X + 8X' Z ) 


dim ^(23) = 2, dimS 4 (23) = 1+ 4 



Here a satisfies a 2 — a + 1 = 0. The rational weight four newform lifts. 



Level 25. rank = 7. 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 + X + 8X 2 ) 


Ella 


2 


T 2 


(1 - 4X)(1 - 8X)(1 + aX + 2X' 2 ) 


EIV 


2 


T 2 


(1 - 2X)(1 - 4X)(1 + PX + 8X 2 ) 


Ellb 


2 


T 2 


(1 -X)(l- 2X)(1 + 7 X + 32X 2 ) 


dim5 2 (25) = 0, dimS 4 (25) = 1 + 1 + 1 



Of the three weight four rational newforms, only one lifts. Here a satisfies a 2 + 1 = 
0, j3 satisfies (3 2 + 49 = 0, and 7 satisfies 7 2 + 16 = 0. 



Level 27. rank =12. 


Ila 


1 


T 2 


(1 -4X)(1-8X)(1 + 2X 2 ) 


lib 


1 


T 2 


(1-X)(1-2X)(1 + 32X 2 ) 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 + 3X + 8X 2 ) 


Ella 


3 


T 2 


(1 - 4X)(1 - 8X)(1 + X)(l + 2X) 


EIV 


3 


T 2 


(1 - 2X)(l - 4X)(1 + X)(1 + 8X) 


Ellb 


3 


T 2 


(1 - X)(l - 2X)(1 + 4X)(1 + BX) 


dimS 2 (27) = 1, dimS 4 (27) = 1 + 1 + 2 



Of the two weight four rational newforms, only one lifts. 
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Level 29. rank = 6. 


Ila 


2 


T 2 
T 3 


(1 - 4X)(1 - 8X)(1 -aX + 2X' 2 ) 
(1 - 9X)(1 - 27X)(1 + aX + 3X 2 ) 


lib 


2 


T 2 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 
(1 - X)(l - 3X)(1 + 9aX + 243X 2 ) 


IV 


2 


n 


(1 - 2X)(1 - 4X)(1 -«X + 8X 2 ) 
(l-3X)(l-9X)(l-(3a + 8)X + 27X 2 ) 


dimS 2 (29) = 2, dimS 4 (29) = 2 + 5 



Here a satisfies a 2 + 2a — 1 = 0. In this example, the weight two and weight four 
newforms that lift are defined over the same quadratic extension of Q. 



Level 31. rank = 6. 


Ila 


2 


T 2 
T 3 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 
(1 - 9X)(1 - 27X)(1 + 2aX + 3X 2 ) 


lib 


2 


T 2 
T 3 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 
(1 - X)(l - 3X)(1 + 18aX + 243X 2 ) 


IV 


2 


T 2 
T 3 


(1 - 2X)(1 - 4X)(1 -0X + 8X 2 ) 
(l-3X)(l-9X)(l + (2/3 + 6)X + 27X 2 ) 


dim5 2 (31) =2, dimS 4 (31) = 2 + 5 



Here a satisfies a 2 — a — 1 = 0, and (5 satisfies j3 2 + 5/9 + 2 = 0. 



Level 33. rank = 10. 


Ila 


1 


T 2 


(1 - 4X)(1 - 8X)(1 - X + 2X 1 ) 


lib 


1 


T 2 


(1 - X) {\ - 2X) {1 -AX + 32X 2 ) 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 + 5X + 8X 2 ) 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 + X + 8X 2 ) 


Ila 


3 




(1 - 4X)(1 - 8X)(1 + 2X + 2X 2 ) 


lib 


3 


T 2 


(1 - X)(l - 2X)(1 + 8X + 32X 2 ) 


dimS 2 (33) = 1, dim5 4 (33) = 1 + 1 + 2 + 2 



The three dimensional eigenspaces are lifts of the weight two newform from level 
11. 



Level 35. rank = 7. 


Ila 


1 


T 2 


(1 -4X)(1-8X)(1 + 2X 2 ) 


lib 


1 


T 2 


(1-X)(1-2X)(1 + 32X 2 ) 


Ila 


2 


T 2 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 


lib 


2 


T 2 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 


IV 


1 


T 2 


(1 - 2X)(1 - 4X)(1 - X + 8X 2 ) 


dim5 2 (35) = 1 + 2, dimS 4 (c 


!5) = 1 + 2 + 3 
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Level 37. rank = 8. 


Ila 


1 


T 2 


(1 - 4X)(1 - 8X)(1 + 2X + 2X 2 ) 


lib 


1 


T 2 


(1 - X)(l - 2X)(1 + 8X + 32X 2 ) 


Ila 


1 


T 2 


(1 -4X)(1-8X)(1 + 2X 2 ) 


lib 


1 


T 2 


(1-X)(1-2X)(1 + 32X 2 ) 


IV 


4 


T 2 


(1 - 2X)(1 - 4X)(1 -«X + 8X 2 ) 


dimS 2 (37) = 1 + 1 


., dimS 4 (37) =4 + 5 



Here a satisfies a A + 6a 3 — a 2 — 16a + 6 = 0. 



Level 39. rank = 10. 


Ila 


1 


T 2 


(1 - 4X)(1 - 8X)(1 - X + 2X 2 ) 


lib 


1 


T 2 


(1 - X)(l - 2X)(1 -AX + 32X 2 ) 


Ila 


2 


T 2 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 


lib 


2 


T 2 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 


IV 


1 


T 2 


(1 -2X)(1-4X)(1 + 8X 2 ) 


IV 


3 


T 2 


(1 - 2X)(1 - 4X)(1 + 5X + 8X 2 ) 


dimS 2 (39) = 1 + 2, dim = 1 + 2 + 2 



Here a satisfies a 2 + 2a — 6 = 0. The three dimensional eigenspaces are lifts of the 
weight four newform from level 13. 



Level 41. rank = 9. 


Ila 


3 


T 2 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 


lib 


3 


T 2 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 


IV 


3 


T 2 


(1 - 2X)(1 - 4X)(1 -PX + 8X 2 ) 


dimS 2 (41) = 3, dim5 4 (41) = 3 + 7 



Here a satisfies a 3 + a 2 — 5a — 1 = 0, and (5 satisfies f3 3 + 3/9 2 — 5/9 — 3 = 0. 



Level 43. rank = 10. 


Ila 


1 


T 2 


(1 - 4X)(1 - 8X)(1 + 2X + 2X 2 ) 


lib 


1 


T 2 


(1 - X)(l - 2X)(1 + 8X + 32X' 2 ) 


Ila 


2 


T 2 


(1 - 4X)(1 - 8X)(1 -aX + 2X 2 ) 


lib 


2 


T 2 


(1 - X)(l - 2X)(1 - AaX + 32X 2 ) 


IV 


4 


T 2 


(1 - 2X)(1 - 4X)(1 -0X + 8X 2 ) 


dim ^(43) = 1 + 2, dimS^ 


13) = 4 + 6 



Here a satisfies a 2 -2 = 0, and (3 satisfies (3 4 + 4(3 3 - 9(3 2 - U(3 + 2 = 0. 
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Level 47. rank = 11. 


Ha 


4 


T 2 


(1 


-4X)(1 


-8X)(1 


-aX + 2X 2 ) 


lib 


4 




(1 


-X){1- 


-2X)(1- 


■ AaX + 32X^) 


Ha 


3 


T 2 


(1 


-4X)(1 


-8X)(1 


-PX + 2X 2 ) 


dimS 2 (47) = 4, dimS 4 (47) = 


= 3 


+ 8 







Here a satisfies a 4 — a 3 — 5a 2 + 5a — 1 = 0, and (3 satisfies /3 3 + 5j3 2 — 2(3 — 12 = 0. 



6. Interpretation of the numerical results 

6.1. The first step, for each Hecke eigenvector (3 G H 5 (T,C) and for each prime I, 
is to write down and factor the Hecke polynomial Pi(X). 

We then see that for the data computed so far, (3 has one of the following Galois 
representations attached. As above, p is a prime not dividing N or any of the Vs we 
are looking at. 

Let k = 2 or 4 and consider the continuous semisimple representation : Gq — > 
GL 2 (Q P ) unramifled outside pX attached to a classical Hecke eigenform / of weight 
k and level N' dividing N. If N is prime, we also assume / has trivial nebentypus. 
Let p = eVfc © t b © e c where (Jfe, a, 6, c) = (2, 0, 2, 3) or (2, 2, 0, 1), or (4, 0, 1, 2). If / 
is an Eisenstein series, then p = xo © Xi e © X2^ 2 © X3£ 3 , where the X; are Dirichlet 
characters of conductor dividing N with values in a finite extension of F, at least two 
of which are trivial. 

Then for any (3 there is some choice of such p which is apparently attached to (3, 
in the sense that the Hecke polynomial at I equals the characteristic polynomial of 
Frobj for all I for which we computed the Hecke eigenvalues. 

Thus it appears that none of our computed classes so far is cuspidal. Therefore we 
should be able to related them to cohomology of the boundary, either geometrically 
or in terms of Eisenstein series. We cannot do this thoroughly, because neither the 
cohomology of the Borel-Serre boundary nor the theory of Eisenstein cohomology has 
been sufficiently worked out for GL 4 / Q. This is not an easy task. We can give the 
following indications. 

From results of Moeglin-Waldspurger [p^| , we don't expect any of our classes to be 
residues of Eisenstein series. In the framework of |17| we can guess that our classes 



lie either in the part of the cohomology indexed by the associate class of parabolic 
subgroups of GL 4 of type (2, 1, 1) or in the part indexed by the Borel subgroup. Here 
the cuspidal data on the GL 2 -factor of the Levi component of the first parabolic comes 
from the appropriate classical cuspform of weight 2 or 4, and we use the appropriate 
power of the determinant on the GLx factors. 

6.2. Geometrically, we make the following comments. Let M be the quotient of 
the symmetric space for SL4(M) by V and let dM be the boundary of its Borel-Serre 
compactification. The covering of dM by its faces gives a spectral sequence for its 
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cohomology. The E 2 page has for its (z, j)-th term fP(Tits building/r, H 3 : (Fiber)). 
An element of that is an assignment: to every face e'(P) of codimension % we assign 
an element of the cohomology in degree j of PflT. These assignments when restricted 
to a common face of codimension i + 1 must add up to 0. Such an assignment gives 
a class in E^ldM). If it persists in the spectral sequence to E^, it will contribute 
to the cohomology H l+J (dM,C). There remains the question as to whether this 
contribution is the restriction of a class in H l+J (M, C). 

Note that if P = LU is a Levi decomposition of P then the spectral sequence of the 
fibration for the cohomology of P fl Y corresponding to this decomposition is known 
to degenerate at E^ q = H p (Tl, H q (U fl T) where is the projection of T to L. 

The classes we have computed so far we expect to be coming in this way from 
dM. From the shape of the apparently associated Galois representations, here is 
what we believe is their origin. We only sketch the constructions, since a detailed 
description would require a thorough investigation of the cohomology of dM for 
arbitrary congruence subgroups of SL 4 (Z). First assume / is a cusp form. 

6.3. The case where a has weight 2: By the Eichler-Shimura theorem, the cuspform 
/ that has a attached shows up as a class a in iJ 1 (A, C), where A is the classi- 
cal Tq(N') C SL/2(Z). First suppose N' = N. Consider the following element of 
i?2' 5 ((9M): O n the standard parabolic subgroup of type (2,2) which is the stabilizer 
of the span of (ei,e2) in 4-space, we put the cohomology class a x 1 on the Levi 
component where we view the trivial coefficients of a as the module H 4 (U fl T). One 
sees that there is a unique class in the appropriate face corresponding to a (3, l)-type 
parabolic subgroup that has the same restriction to the type (2, 1, 1) face they have 
in common and restricts to on the other faces. Hence these two glue up to give 
a class in the boundary. The same construction with the transposed parabolic sub- 
group also gives a class, and these seem to account for all our /3's falling under this 
case. If N' ^ N, we choose a (2, 2)- parabolic subgroup P such that T L has level N'. 
Then we imitate the construction above. 

6.4. The case where a has weight 4: Here not every class we construct on the 
boundary seems to lift to if 5 (M), but only some of them. We don't know the reason 
for this. The construction here creates a class in El' 4 (dM), and if this correctly 
describes what we have computed, our computed classes of this type would be ghost 
classes. That is, the corresponding cohomology class in H 5 (M) restricts nontrivially 
to the boundary of M, but it restricts to on each face of the boundary, since it is 
coming from a class in E^^dM) with % > 0. 

For this construction, one first chooses a parabolic subgroup P of type (2,1,1). 
Note that H 3 ((UnT), C) contains a T^-submodule isomorphic to V 2 , the homogeneous 
complex polynomials of degree 2 on 2-space, after we identify T L with a subgroup of 
GL 2 (Z)xGL 1 (Z)xGL 1 (Z). By the Eichler-Shimura theorem, the weight 4 cuspform / 
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that has a attached shows up in H 1 (A, V2). Thus we can view it as in H 1 (Tl, H 3 (U n 

H). 

We can do this on three P's which are not conjugate to each other in such a way 
that they, together with certain classes on (3, l)-type parabolic subgroups, all glue 
together to give a class in E^(dM). The details are left for the reader. 

6.5. Finally, the case where a is an Eisenstein series is harder to understand. We 
haven't worked out exactly how the gluing process goes in this case, so we're not sure 
what stratum of the spectral sequence is occupied by the corresponding boundary 
classes. 

We note that in every case, p restricted to an inertia subgroup of Gq at p has the 
form 1 © e © e 2 © e 3 , which is consistent with the conjecture of Ash-Sinnott [10 



since 



the coefficient module of the cohomology classes we consider is the trivial module. 
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